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We calculate the four-loop /3-function for the generalised Wess-Zumino model. We 
use the result and Pade-Borel summation to discuss the domain of attraction of the quasi- 
infra-red fixed point of the top-quark Yukawa coupling in the supersymmetric standard 
model, and argue that the domain is in fact substantial. 
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In this note we present the four-loop /^-function for a many-field Wess-Zumino model 
with arbitrary cubic interactions. (This calculation is a straightforward generalisation of 
the existing one of Avdeev et al|l[ for the single field and coupling case). We use the result 
to write down the 0(\®) terms in the f3- function for the top-quark Yukawa coupling in the 
supersymmetric standard model. We then use this result to discuss issues regarding the 
limitations of perturbation theory in the standard running analysis whereby the low energy 
supersymmetric standard model is matched onto a much reduced (in terms of the number 
of free parameters) theory at high energies. In such analyses, putative supersymmetric 
spectra are often presented as scatter-plots, for ranges of certain of the free parameters, 
subject to certain cuts. One such cut, frequently applied, is that the Yukawa couplings 
remain perturbative throughout: in for example Ref. one finds that \t,b,r < 3.5 has been 
enforced. What is the basis and reliability of such cuts, given the (presumably) asymptotic 
nature of the perturbation expansion? This question becomes particularly interesting in 
the context of the possible quasi-infra-red fixed point (QIRFP) behaviour of At 

We begin with the generalised Wess-Zumino model, defined by the superpotential: 

W = \Y ijk ^^ k . (1) 

where $i is a multiplet of chiral superfields. The /3-functions for the Yukawa couplings 
/3y fc are given by 

0»k = yj>(« 7 k) p = y^ k p + ( k ^l) + (k~ j), (2) 

where 7 is the anomalous dimension for <E>. 

The result for 7 through three loops is as follows 0: 

16tt 2 7 (1) = P (3a) 
(16tt 2 ) 2 7 (2) = -Si (36) 
(16tt 2 ) 3 7 (3) = f C(3)M + 2Y*S 4 Y - \S 7 - S 8 (3c) 

where our notation follows that of Ref. , except that here we have no gauge coupling: 

pi . = i Y ikl Y jkl (4a) 
S\j = Y* mn pP m Y Jpn (46) 
{Y*S A Y) i j = Y imn S 4 P m Y wn . (4c) 

qi . — y imn PP PI V- (Arl\ 
u 7j 1 1 m 1 n ± jpq \ "'J 

Sl 3 =Y* mn (P 2 y m Y 3pn (4e) 
M*, = Y M Y kmn Y lrs YP mr Y^ s Y Jpq . (4/) 
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Our result for 7^) is as follows: 

(167r 2 )y 4 ) = lY* km S 7 l k Y jlm + (f - 2C(3))Y* km S 8 l k Y jlrn + ±Y ikm {PS A + S 4 P) l k Y jlm 

- 5Y^Y kmp S 4 n m Y ln PY jlm - (f C(3) + ^W lkm M l k Y jlm 
+ (2((3) - 1) [Y ikm (P 3 ) l k Y jlm +Y* km (P 2 ) l k P n rn Y jln ] 

+ lY ikm P l k (S^) n m Y jln - 10((5)Y M Y kmn Y lsv Y nsu Y rtu Yr mr Y^Y Jpq 

- (6C(3) - 3C(^MY tkt P l t Y kmn Y lrs YP mr Y^ s Y Jpq 

1 lyiHy y ypmrynstpg y iV ikl V, p t V, -ypmr-yqns-y l 
~r 2 1 kmn 1 lrs 1 1 1 t 1 jpq > 1 1 kmt 1 n 1 lrs 1 1 3V1\' 

(5) 

We have explicitly calculated the requisite Feynman diagrams for the most part. (This 
was in any case necessary since in Ref. results are given for subsets of, not individual 
graphs.) The results for several of the Feynman integrals are given in Ref. 0. In two 
places, however, we have relied on previous authority: specifically, in Table 1 of Ref. 
we have used their total for the set of graphs 4. 1 to finesse the calculation of one particular 
graph; and we have used the result of Ref. for the Feynman integral that arises in the 
calculation of the graph 4.5. 

For the supersymmetric standard model superpotential, which is (retaining only At) 

W = X t H 2 Qt (6) 
where Q = , we find, using (3\ t = (jh 2 + 1q + 7t) X t M : 

16n 2 p[] } = 6A t 3 (7a) 

(IGtt 2 ) 2 /^ = -22X 5 t (76) 

(167r 2 ) 3 /3g ) = [102 + 36C(3)]Aj (7c) 

and finally from Eq. (|5|): 

(16n 2 ) 4 p[^ = - [678 + 696C(3) - 216C(4) + 1440C(5)] X 9 t . (8) 

A specific aspect of the running analysis that has become popular in recent years 
is the possibility that At may exhibit QIRFP behaviour. The main attraction of this 
philosophy is the idea that a large range of input couplings at high energies may produce 
the same value of At at Mz- It is clearly interesting to inquire as to the extent to which 
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this range is limited by the requirement of perturbative believability. Let us review the 
QIRFP paradigm: the one-loop equations governing the evolution of the At and the gauge 
couplings are 

d(Xi 



d Vt /a 16 

-jj: = yt{6yt ~ -fa3 



3«2 — ii a i)i an d 



dt 



biOii 



(9) 



where y t = ^A t , 



on 



«1,2,3, h 
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1 —3), and t = ^-ln/i. We have made the 
approximation that y t » y T , yj,, which involves the assumption that it is not the case that 
tan/3 >> 1; it is straightforward to consider the more general case but we omit it here for 
simplicity. From Eq. (|9|) we obtain: 



3/t(M 



Xj 



yt{M 2 



/(0) - 6Fy t (M z ) 



where 



f(r)dr, T 



1 , M x 
— In 

2tt M z 



and 



fir) 



r« 3 (T)" 


16 
3i> 3 


~a 2 (TY 


3 


'ai(T)" 


13 
15b 1 






. a 2(r) _ 




.«i( T ) . 





(10) 



(12) 



We will take Mx to be the unification scale, requiring that a 2 (T) = a\(T). Of course 
it is a straightforward matter to integrate the differential equations numerically, but the 
partial analytic solution above is nevertheless useful in order to see what is going on. From 
Eq. (Jl(]) we see that y t suffers a Landau pole unless 



y t (M z ) < /(0)/(6F). 
Another nice way of thinking about this result is as follows. Suppose 

yt(Mx) » 1/(6F). 
Then it is easy to see from Eq. (|T0| ) that in this limit 

y t (M z )=y t = f(0)/(6F) 



(13) 



(14) 



(15) 



which is independent of j/t(Mi)! Of course y t is the value of yt(Mz) such that the Landau 
pole occurs at Mx, and hence represents an upper limit; but it is more productive to think 
of yt(Mx) as the input, as follows. Depending on the value of F, there may be a wide range 
of values of yt(Mx) which all lead to the same value of yt(Mz)- The corresponding value, 
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Vt(Mz) = /(0)/(6F), is called a quasi-infra-red fixed point of the evolution equations. The 
prefix quasi- is occasioned by the fact that f/(6F) is of course a function of T. This 
behaviour is illustrated in Fig. 1. 
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Fig.l: Plot of y t against log 10 lc f eV for various values of yt(Mx) 

If we take the limit T — > oo, then we approach genuine fixed point behaviour H. It 
is easy to see this in the approximation a± = a.2 = 0, for which the equations for -j^ and 
^f- exhibit an infra-red fixed point such that 



a 3 



-3 + 



16 



6 



0.39. 



For a 3 (M z ) = 0.1 this gives y t (M z ) = 0.039. Substituting in Eq. (|TT|) and (0) (with 
a x {Mz) = 0.0167, a 2 {M z ) = 0.0320) one finds /(0) = 9.98 and F = 17.88 so that 
y t = 0.093. Thus QIRFP and IRFP behaviour are quite distinct; it is the former which 
is relevant for the supersymmetric standard model. (This was first pointed out for the 
standard model in Ref. 

We now see that since 1/(6F) w 0.01 we have that there is a wide range for yt(Mx) 
such that 

1/(6F) « y t (M x ) < 1 (16) 

where the upper limit is a naive constraint for perturbative believability. (Note that the 
constraint A < 3.5 used in Ref. would correspond to yt < 0.975). We want to address 
this latter restriction; does our four-loop calculation afford any insight into it? In the 
literature, cuts of the type At < 3.5 alluded to above are sometimes motivated by requiring 
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f3^ < xf3^ l \ where is the appropriate L-loop /3-function, and x is some convincing 
fraction: | for example. This leads, using Eq. ([F[) to yt < 0.86. A naive extension of this 
approach would obviously suggest a drastic curtailment of the acceptable range of yt(Mx), 
since for instance (3^ < gives yt < 0.16 0. This is illustrated in Fig. 2, where we 

plot y t {Mx) against y t (M z ) for L = 1, • • -4. 
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Fig. 2: Plot of yt(Mx) against y t (Mz) ■ The solid, dashed, dotted and 
dash-dotted lines correspond to one, two, three and four-loop (3-functions 
respectively. 

We wish to argue that the (presumably) asymptotic nature of the perturbation series 
for the /3-functions means that in fact the actual domain of attraction of the QIRFP is 
more accurately represented by the one-loop than by the four-loop approximation. 

A striking feature of Eq. (|[) is the broad similarity to the corresponding result for 
0(n) 4> . Most importantly, note the characteristic alternating sign behaviour, suggesting 
the possibility of Borel summability. (For a review on the problem of resummation of 
perturbation theory see Ref. p^.) We have not found in the literature any discussion of the 
large order behaviour of the Wess-Zumino model; but for the supersymmetric anharmonic 



oscillator, it has been noted [14] that while the supersymmetric case does represent a 



Such naive constraints may be even more restrictive in other sectors of the theory: in 



particular the soft supersymmetry breaking sector. By demanding that (3 \ < -r(3 2 • 

m Q m Q 

yt < 0.32, for any theory with the commonly assumed universal form of the soft supersymmetry 



we get [|1( 
symmetr 

breaking parameters (specific models might originate smaller upper bounds; the "P = class 



of models [11], for example, gives yt < 0.21). 
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bifurcation point with respect to some behaviour, nevertheless the characteristic factorial 
divergence L\ persists. It is interesting to note, however, that the exact form for the gauge 
/9-function in an N 



1, d = 4 theory without chiral fields, [15j [|16 

-3C(G) 



ft, 



9 
16tt s 



1 - 2C(G)g 2 (16ir 2 ) 



(17) 



clearly has a finite radius of convergence. Now this result may not hold in DRED (di- 
mensional reduction with minimal subtraction); it was shown in Ref. [17] that the scheme 



in which Eq. (17) is valid differs from DRED. In any case, for the Wess-Zumino model, 
there is certainly no indication from Eq. (|5|) that the Yukawa 8- function is other than 
an asymptotic Lies! Now the existence of _>„n singnlanties H implies that 
in asymptotically un-free theories (such as the Wess-Zumino model) amplitudes are not 
Borel-summable in general; but specific renormalisation group functions may be so. It 
seemed to us worthwhile exploring the consequences of naive Pade-Borel (PB) summation 
of the yt-dependence of (3 yt . 

PB summation proceeds as follows (see for example Ref. [M). Given a series 



f( x ) = 5Z ar 



x 



n=0 



one defines 



71 = 



and calculates [N,M] Pade approximants to B(x), Bnm- 
of f(x) is given by 



(18) 
(19) 

Then the PB-summed version 



F N m(x) 



e t B NM (xt)dt. 



(20) 



Essentially this construction amounts to a guess of the coefficients of powers of x beyond 
those originally calculated, incorporating (for x L ) a factor L\. 

We have calculated the [1, 1] and the [2, 2] Pades for the t- Yukawa /^-function, i.e. the 
series 

4^22^ 2 + (^) 2 [102 + 36C(3)]y t 3 

(21) 

(_L)3 [ 678 + 696C(3 ) _ 216C(4) + 1440C(5)] yf + 



^ = y t {0 + 6y t 



2 Note that if this was true the same would hold for the exact result of Ref. []l6|| for (3 g , when 
chiral fields are present. 
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We have written the series in this form in order that the 0(y^) term in the curly brackets 
should correspond to the Lth order contribution in perturbation theory. The L\ factor 
produced by the PB process for the 0(y^) term in the series then correctly mimics the 
expected LI growth of the L-loop perturbation theory contributions. The [1, 1] PB approx- 
imation for the series 



dyt 
dt 



y t (0 + ay t -by\ 



is easily calculated as 



dyt 
dt 



Vt 



2a 2 



4a 3 eT ^ E / 2a 
b 2 yt 1 \by t 



(22) 



(23) 



where Ei(x) = ^-j-dt is the exponential integral function. The [2, 2] PB approximation 
for the series 



dyt 
dt 



y t (0 + ay t - by 2 + cyf - dy\ 



is found to be 

dy t 



dt 
where 



= Vt 



rfr 2 (ar 1 -a) _n_ / rA r x r\(ar 2 -d) _n / r 2 

vrir 2 + — t e w EW + — r e * Ei 

(r 1 -r 2 )yt \ yt J [r 2 -r 1 )y t \ y t 



a 2 d + 3b s - Aabc 
4ac — 6b 2 



(24) 



, (25) 



(26) 



and r% and r 2 are the roots of 



(3bd - 4c 2 )r 2 + (Gad - 12bc)r + 12(2ac - 3b 2 ) = 0. 



(27) 



In each case the PB approximation involves the exponential integral function; this is a 
generic feature. Note that with a . . .d as given in Eq. ( pl|) , the roots of Eq. ( p7|) are 
both negative and real. It follows that in Eq. ( |25|) , no singularity is encountered in the 
exponential integrals; the same is true of Eq. (^3|). In the [1,2] PB, however, the roots of 
the analogous quadratic have opposite signs. Now such a pole may represent a physical 
renormalon singularity; but since we do not encounter it in the [2, 2] Pade, and because 
we have no specific knowledge of the asymptotic behaviour for this theory, we choose to 
ignore this case. 

We have investigated the running of yt and cti,2,3 by adding the perturbative contribu- 
tions involving the gauge couplings to Eq. (0) and evolving the gauge couplings using the 
perturbative /3-functions. This seems justified on the grounds that cti 5 2,3 remain small over 
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the range of interest. The results are displayed in Fig. 3, where we plot y t (M x ) against 
Ut{M z ) for our [1,1] PB and [2,2] PB approximations. For purposes of comparison, we 
also give the results obtained using the one-loop and four-loop /3-functions. (Since we only 
know explicitly the pure Yukawa contribution at the four-loop level, we must omit the 
gauge contributions at this level; but in the region of interest we know that the Yukawa 
contributions are dominant.) 

0.096 
0.094 
0.092 
0.09 
0.088 

N 

1 0.086 
0.084 
0.082 
0.08 
0.078 
0.076 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

yt at MX 

Fig. 3: Plot ofyt(Mz) against yt(Mx)- The solid and dotted lines cor- 
respond to the one and four-loop perturbative (3-functions, and the dot- 
dashed and dashed lines to the [1, 1] and [2, 2] PBs respectively. 

For small y t: the [1, 1] and [2,2] PB results should approach the two-loop and four- 
loop perturbative results respectively; this is clearly seen in the [2,2] case. For large yt, 
the asymptotic nature of the perturbation series implies that lower orders in perturbation 
theory should be more accurate than higher orders; accordingly, we see that as yt increases, 
the PB evolution starts to resemble the one-loop perturbative behaviour. Furthermore, 
the [2,2] PB is closer to the one-loop result than the [1, 1]. We have also repeated this 
exercise for the [2, 1] PB; the results are almost indistinguishable from those for the [2, 2] 
Pade and indicate that the successive PBs may be converging quite rapidly in the region 
of interest. 

The one-loop perturbative evolution displays the QIRFP behaviour as explained in our 
earlier analysis; we see, for instance that values of y t (M x ) in the range 0.2 < y t (M x ) < 1 
lead to values of y t (M z ) in the range 0.089 < y t (M z ) < 0.092. The upper limit on 
y t {Mx), as discussed earlier, is a rough constraint for perturbative believability, and if 
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one wished to include higher-loop perturbative contributions there is a prima facie case 
for restricting the allowed range of y t (Mx) still further. The PB results show the same 
QIRFP behaviour; however, we believe they are reliable up to a larger value of y t (M x )- 
Of course, there must still be an upper limit on y t beyond which we can no longer trust 
the PB approximation. Some indication of where this upper limit is can be obtained by a 
comparison with alternative means of implementing the PB programme; clearly one can no 
longer trust the PB results in regions where different approaches give qualitatively different 
behaviour. For instance, one could write 



dy t 2 



1 '^fvt + |[102 + 36C(3)K 2 + . . . (28) 



and then perform the PB procedure on the series 1 + ay t + byf + cyf + — Note that 
the L-loop contribution to the perturbation expansion corresponds to the 0(yf J ~ 1 ) term 
in this series. This appears at first to be a drawback, as the PB process will not now 
exactly match the LI growth of the L-loop perturbation theory contribution. However we 
might well expect behaviour of the L-loop contribution like L!L a , and our uncertainty as 
regards the value of a means that we cannot really distinguish between growth like L! or 
like (L- 1)!. 

If one evolves y t alone (setting 0^2,3 = Ub = Ut = 0), using the [1,1] or [2,1] PB 
obtained in this way, one finds a fixed point at large y t , for y t ~ 100 and y t ~ 8 respectively, 
in contrast with the Landau pole behaviour of the previous PBs. However, in the presence 
of the other couplings, the behaviour is somewhat modified. In fact, the evolution is 
practically indistinguishable from that for the [2, 2] PB displayed in Fig. 3 over the range 
shown. Moreover, for larger y t (Mx), the fixed point is obscured by the development of 
a Landau pole in y^ and by the increasing size of yf as terms which have not been Borel 
summed. Nevertheless, the last vestige of a fixed point is discernible in the [2, 1] PB 
evolution at yt ~ 6. This may perhaps be taken as a signal that we should not trust any 
of our PBs beyond this point. Nevertheless we have considerably extended the domain of 
attraction of the QIRFP as compared to the perturbative case beyond y t (Mx) = 1, and 
perhaps optimistically as far as yt(Mx) ~ 5; we now see that values of yt(Mx) in the 
range 0.2 < y t (M x ) < 5 lead to values of y t (M z ) in the range 0.089 < y t (M z ) < 0.094. 
Note that the fixed point for the [2, 2] PB is roughly 2% higher than the effective upper 
limit on y t (Mz) which applies in the one-loop perturbative case; for a fixed value of nit, 
this leads to a reduction in tan/? of around 4.5%. 
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In conclusion, we have demonstrated by means of Pade-Borel summation that the 
domain of attraction of the quasi-infra-red fixed point in the supersymmetric standard 
model is in fact large. There has been some recent speculation |2(| on the possible role 
of QIRFP behaviour in the theory above Mx in assuring universality of the soft-breaking 
parameters. Because the energy range between Mpi anc k and Mx is much smaller than 
that between Mx and Mz, rapid evolution of the couplings is essential if a QIRFP is to 
be approached in this region. Since the approach to the QIRFP appears to be quicker for 
larger initial couplings, it might be of interest to use PB summation to explore the region 
of larger coupling with more confidence. 
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